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Abstract. We present a review of vector field models of inflation and, in particular, 
of the statistical anisotropy and non-Gaussianity predictions of models with SU (2) 
vector multiplets. Non-Abelian gauge groups introduce a richer amount of predictions 
compared to the Abclian ones, mostly because of the presence of vector fields self- 
interactions. Primordial vector fields can violate isotropy leaving their imprint in the 
comoving curvature fluctuations C at late times. We provide the analytic expressions 
of the correlation functions of £ up to fourth order and an analysis of their ampli- 
tudes and shapes. The statistical anisotropy signatures expected in these models are 
important and, potentially, the anisotropic contributions to the bispectrum and the 
trispectrum can overcome the isotropic parts. 
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1. Introduction 



In the standard cosmological model, at very early times the Universe undergoes a quasi 
de Sitter exponential expansion driven by a scalar field, the inflaton, with an almost 
flat potential. The quantum fluctuations of this field are thought to be at the origin of 
both the Large Scale Structures and the Cosmic Microwave Background (CMB) fluc- 
tuations that we are able to observe at the present epoch [1]. CMB measurements 
indicate that the primordial density fluctuations are of order 10 -5 , have an almost 
scale-invariant power spectrum and are fairly consistent with Gaussianity and statisti- 
cal isotropy [2, 3, 4, 5, 6]. All of these features find a convincing explanation within 
the inflationary paradigm. Nevertheless, deviations from the basic single-(scalar)field 
slow-roll model of inflation are allowed by the experimental data. On one hand, it is 
then important to search for observational signatures that can help discriminate among 
all the possible scenarios; on the other hand, it is important to understand what the 
theoretical predictions are in this respect for the different models. 

Non-Gaussianity and statistical anisotropy are two powerful signatures. A random 
field is defined "Gaussian" if it is entirely described by its two-point function, higher 
order connected correlators being equal to zero. Primordial non-Gaussianity [7, 8] 
is theoretically predicted by inflation: it arises from the interactions of the infla- 
ton with gravity and from self-interactions. However, it is observably too small in 
the single-field slow-roll scenario [9, 10, 11]. Alternatives to the latter have been 
proposed that predict higher levels of non-Gaussianity such as multifield scenarios 
[12, 13, 14, 15, 16, 17, 18], curvaton models [19, 20, 21, 22, 23, 24] and models with 
non-canonical Lagrangians [25, 26, 27, 28, 29]. Many efforts have been directed to the 
study of higher order (three and four-point) cosmological correlators in these models 
[30, 33, 34, 31, 28, 11, 32, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45] and towards im- 
proving the prediction for the two-point function, through quantum loop calculations 
[46, 47, 48, 10, 49, 50, 51, 52, 53]. From WMAP, the bounds on the bispectrum amplitude 
are given by -4 < f§ c L < 80 [54] and by -125 < < 435 [55] at 95% CL, respectively 
in the local and in the equilateral configurations. For the trispectrum, WMAP provides 
—5.6 x 10 5 < gNL < 6.4 x 10 5 [56] (gNL is the "local" trispectrum amplitude from cubic 
contributions), whereas from Large-Scale-Structures data —3.5 x 10 5 < g^L < 8.2 x 10 5 
[57], at 95% CL. Planck [58] and future experiments are expected to set further bounds 
on primordial non-Gaussianity. 

Statistical isotropy has always been considered one of the key features of the CMB 
fluctuations. The appearance of some "anomalies" [59, 60, 61] in the observations 
though, after numerous and careful data analysis, suggests a possible a breaking of 
this symmetry that might have occurred at some point of the Universe history, possibly 
at very early times. This encouraged a series of attempts to model this event, preferably 
by incorporating it in theories of inflation. Let us shortly describe the above mentioned 
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"anomalies" . First of all, the large scale CMB quadrupole appears to be "too low" and 
the octupole "too planar"; in addition to that, there seems to exist a preferred direction 
along which quadrupole and octupole are aligned [62, 63, 59, 64, 65]. Also, a "cold spot" , 
i.e. a region of suppressed power, has been observed in the southern Galactic sky [60, 66]. 
Finally, an indication of asymmetry in the large-scale power spectrum and in higher- 
order correlation functions between the northern and the southern ecliptic hemispheres 
was found [67, 61, 68]. Possible explanations for these anomalies have been suggested 
such as improper foreground subtraction, WMAP systematics, statistical flukes; the 
possibilities of topological or cosmological origins for them have been proposed as well. 
Moreover, considering a power spectrum anisotropy due to the existence of a preferred 
spatial direction h and parametrized by a function g(k) as 



the five-year WMAP temperature data have been analyzed in order to find out what 
the magnitude and orientation of such an anisotropy could be. The magnitude has been 
found to be g = 0.29 ± 0.031 and the orientation aligned nearly along the ecliptic poles 
[69]. Similar results have been found in [70], where it is pointed out that the origin of 
such a signal is compatible with beam asymmetries (uncorrected in the maps) which 
should therefore be investigated before we can find out what the actual limits on the 
primordial g are. 

Several fairly recent works have taken the direction of analysing the consequences, in 
terms of dynamics of the Universe and of cosmological fluctuations, of an anisotropic 
pre-inflationary or inflationary era. A cosmic no-hair conjecture exists according to 
which the presence of a cosmological constant at early times is expected to dilute any 
form of initial anisotropy [71]. This conjecture has been proven to be true for many (all 
Bianchi type cosmologies except for the the Bianchi type- IX, for which some restrictions 
are needed to ensure the applicability of the theorem), but not all other kinds of metrics 
and counterexamples exist in the literature [72, 73, 74]. Moreover, even in the event 
isotropization should occur, there is a chance that signatures from anisotropic inflation 
or from an anisotropic pre-inflationary era might still be visible today [75, 76, 77, 78]. 
In the same contest of searching for models of the early Universe that might produce 
some anisotropy signatures at late time, new theories have been proposed such as spinor 
models [79, 80, 81, 82], higher p-forms [83, 84, 85, 86, 87, 88] and primordial vector field 
models (see Section 2 for a quick review). 

We are going to focus on statistical anisotropy and non-Gaussianity predictions of pri- 
mordial vector field models. As mentioned above, there are great expectations that 
Planck and new experiments will, among other things, shed more light on the level of 
non-Gaussianity of the CMB fluctuations and on the nature of the unexpected anisotropy 
features we mentioned (see, e.g., [89]). Models that combine both types of predictions 
could be more easily testable and, from non-Gaussianity measurement, more stringent 
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statistical anisotropy predictions could be produced or viceversa. 

Within vector field models, higher order correlators had been computed in [90, 91, 92, 
93, 94] and, more recently, in [95, 96] for U(l) vector fields. We considered SU(2) vec- 
tor field models in [97, 98]. Non-Abelian theories offer a richer amount of predictions 
compared to the Abelian case. Indeed, non-Abelian self interactions provide extra con- 
tributions to the bispectrum and trispectrum of curvature fluctuations that are naturally 
absent in the Abelian case. We verified that these extra contributions can be equally 
important in a large subset of the parameter space of the theory and, in some case, can 
even become the dominant ones. 

This paper is structured as follows: in Sec. 2 we review some vector field models of 
inflation; in Sec. 3 we present the SU (2) model; in Sec. 4 we provide the results for the 
two, three and four-point functions of the curvature fluctuations; in Sec. 5 we present 
the non-Gaussianity amplitudes for the bispectrum and for the trispectrum; in Sec. 6 
we show and discuss their shapes; finally in Sec. 7 we draw our conclusions. 

2. Inflation and primordial vector fields 

The attempt to explain some of the CMB "anomalous" features as the indication of a 
break of statistical isotropy is the main reason behind ours and many of the existing 
inflationary models populated by vector fields, but not the only one. The first one of 
these models [99] was formulated with the goal of producing inflation by the action of 
vector fields, without having to invoke the existence of a scalar field. The same motiva- 
tions inspired the works that followed [100, 101, 102]. Lately, models where primordial 
vector fields can leave an imprint on the CMB have been formulated as an alternative to 
the basic inflationary scenario, in the search for interesting non-Gaussianity predictions 
[90, 91, 92, 93, 94, 97, 98, 95, 96]. Finally, vector fields models of dark energy have 
been proposed [103, 104, 105, 106, 107, 108]. All this appears to us as a rich bag of 
motivations for investigating these scenarios. 

Before we quickly sketch some of them and list the results so far achieved in this di- 
rection, it is important to briefly indicate and explain the main issues and difficulties 
that these models have been facing. We will also shortly discuss the mechanisms of 
production of the curvature fluctuations in these models. 

Building a model where primordial vector fields can drive inflation and/or produce the 
observed spectrum of large scale fluctuations requires a more complex Lagrangian than 
the basic gauge invariant C vector = —(y/^g/^F^F^. In fact, for a conformally in- 
variant theory as the one described by £ vector, vector fields fluctuations are not excited 
on superhorizon scales. It is then necessary to modify the Lagrangian. For some of 
the existing models, these modifications have been done to the expense of destabilizing 
the theory, by "switching on" unphysical degrees of freedom. This was pointed out in 
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[109, 110, 111], where a large variety of vector field models was analyzed in which lon- 
gitudinal polarization modes exist that are endowed with negative squared masses (the 
"wrong" signs of the masses are imposed for the theory to satisfy the constraints that 
allow a suitable background evolution). It turnes out that, in a range of interest of the 
theory, these fields acquire negative total energy, i.e. behave like "ghosts", the presence 
of which is known to be responsible for an unstable vacuum. A related problem for some 
of these theories is represented by the existence of instabilities affecting the equations 
of motion of the ghost fields [109, 110, 111]. 

In the remaing part of this section, we are going to present some of these models to- 
gether with some recent attempts to overcome their limits. 

In all of the models we will consider, primordial vector fields fluctuations end up either 
being entirely responsible for or only partially contributing to the curvature fluctuations 
at late times. This can happen through different mechanisms. If the vector fields affects 
the universe expansion during inflation, its contribution Qa to the total ( can be derived 
from combining the definition of the number of e-foldings (N = J Hdt) with the Einstein 
equation (H 2 = (8iiG/3)(p ( f > + pa), Pa being the energy density of the vector field and 
p^ the inflaton energy density) and using the <5N expansion of the curvature fluctuation 
in terms of both the inflaton and the vector fields fluctuations (see Sec. 4). To lowest 
order we have [92] 

Ca = t^tSA, (2) 
2m P 

where a single vector field has been taken into account for simplicity (mp is the reduced 
Planck mass, A is the background value of the field and 5 A its perturbation). When 
calculating the amplitude of non-Gaussianity in Sec. 5, we will refer to this case as 
"vector inflation" for simplicity. 

A different fluctuation production process is the curvaton mechanism which was initially 
formulated for scalar theories but it is also applicable to vectors [112, 113]. Specifically, 
inflation is driven by a scalar field, whereas the curvaton field(s) (now played by the 
vectors), has a very small (compared to the Hubble rate) mass during inflation. Towards 
the end of the inflationary epoch, the Hubble rate value starts decreasing until it equates 
the vector mass; when this eventually happens, the curvaton begins to oscillate and it 
will then dissipate its energy into radiation. The curvaton becomes responsible for a 
fraction of the total curvature fluctuation that is proportional to a parameter (r) related 
to the ratio between the curvaton energy density and the total energy density of the 
universe at the epoch of the curvaton decay [92] 
, r bp A 

Ca - -— , (3) 
<j Pa 

where r = 3pa/{3pa + 4p</,). Anisotropy bounds on the power spectrum favour small 
values of r. 

From Eqs. (2) and (3) we can see that, dependending on which one of these two mecha- 
nisms of production of the curvature fluctuations is considered, different coefficients will 
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result in the SN expansion (see Eq. (22)). 



In this section we will describe both models where inflation is intended to be vector field 
driven and those models in which, instead, the role of the inflaton is played by a scalar 
field, whereas the energy of the vector is kept subdominant in the total energy density 
of the universe during the entire inflationary phase. 

2.1. Self-coupled vector field models 

A pioneer work on vector field driven inflation was formulated by L. H. Ford [99], who 
considered a single self-coupled field with a Lagrangian 

Uector = —F^F^ + Vty) (4) 

where = d fl B u — d u B fl and the potential V is a function of ip = B a B a . Different sce- 
narios of expansion are analyzed by the author for different functions V. The universe 
expands anisotropically at the end of the inflationary era and this anisotropy either 
survives until late times or is damped out depending on the shape and the location of 
the minima of the potential. 

The study of perturbations in a similar model was proposed by Dimopoulos in [112] 
where he showed that for a Lagrangian 

Uector = --F^ V + l -m 2 B^ (5) 

and for m 2 ~ —2H 2 , the transverse mode of the vector field is governed by the same 
equation of motion as a light scalar field in a de Sitter stage. A suitable superhorizon 
power spectrum of fluctuations could therefore arise. In order to prevent production of 
large scale anisotropy, in this model the vector field plays the role of the curvaton while 
inflation is driven by a scalar field. 

2.2. Vector-field coupled to gravity 

The Lagrangian in Eq. (5) may be also intended, at least during inflation, as including 
a non-minimal coupling of the vector field to gravity; indeed the mass term can be 
rewritten as 

Uector D \ (m\ + iB) B^B" (6) 

where, for the whole duration of the inflationary era, the bare mass m is assumed to 
be much smaller than the Hubble rate and the Ricci scalar R = —6 - + f-1 can be 
approximated as R ~ —12H 2 . For the specific value £ = 1/6, Eq. (5) is retrieved. 

For the Lagrangian just presented, Golovnev et al [100] proved that the problem of 
excessive anisotropy production in the case where inflation is driven by vector fields can 
be avoided if either a triplet of mutually orthogonal or a large number N of randomly 
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oriented vector fields is considered. 



The Lagrangian (6) with £ = 1/6 was also employed in [113], where inflation is scalar- 
field-driven and a primordial vector field affects large-scale curvature fluctuations and, 
similarly, in [114], which includes a study of the backreaction of the vector field on the 
dynamics of expansion, by introducing a Bianchi type-I metric. 

2.3. Ackerman-Carroll-Wise (ACW) model 

A model was proposed in [115] where Lagrange multipliers (A) are employed to determine 
a fixed norm primordial vector field B^B^ = m 2 



where p\ is a vacuum energy. The expansion rate in this scenario is anisotropic: if we 
orient the x-axis of the spatial frame along the direction determined by the vector field, 
we find two different Hubble rates: along the x-direction it is equal to 



and it is given by H a = (1 + cfi 2 )Hb along the orthogonal directions; \i = m/mp, P is 
a polynomial function of [i and c is a parameter appearing in the kinetic part of the 
Lagrangian that we omitted in (7) (see [115] for its complete expression). As expected, 
an isotropic expansion is recovered if the vev of the vector field is set to zero. 

2-4- Models with varying gauge coupling 

Most of the models mentioned so far successfully solve the problem of attaining a slow- 
roll regime for the vector-fields without imposing too many restrictions on the param- 
eters of the theory and of avoiding excessive production of anisotropy at late times. 
None of them though escapes those instabilities related to the negative energy of the 
longitudinal modes (although a study of the instabilities for fixed-norm field models was 
done in [116] where some stable cases with non-canonical kinetic terms were found). As 
discussed in [109, 110, 111], in the self-coupled model a ghost appears at small (com- 
pared to the horizon) wavelengths; in the non-minimally coupled and in the fixed-norm 
cases instead the instability concerns the region around horizon crossing. 

Models with varying gauge coupling can overcome the problem of instabilities and have 
recently attracted quite some attention. In [90], the authors consider a model of hybrid 
inflation [117, 118, 119, 120] with the introduction of a massless vector field 



where <p is the inflaton and x is the so-called "waterfall" field. The potential V is chosen 
in such a way as to preserve gauge invariance; this way the longitudinal mode disappears 




(7) 



L D l - (d^<f> + dpX&x) ~ \f\<P)F^ v + V(<f>, X , B,) 



(9) 
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and instabilities are avoided. 

Similarly, Kanno et al [121] consider a vector field Lagrangian of the type 

Kector = ~f 2 (<f>)F^F^, (10) 

but in a basic scalar field driven inflation model. Very recently, in [122] the linear 

perturbations in these kind of models have been investigated. 

Finally, in [92, 96] varying mass vector field models have been introduced 

L vect or = -\f{<p)F^F^ + ^m 2 B^, (11) 

where / ~ a a and m ~ a (a is the scale factor and a is a numerical coefficient). The 
special cases a = 1 and a = —2 are of special interest. In fact, introducing the fields 
An and Ap, related to one another by A^ = fB^ = aA^ (A^ and A^ are respectively 
the comoving and the physical vectors), it is possible to verify that the physical gauge 
fields are governed by the same equations of motion as a light scalar field in a de Sitter 
background. Vector fields in this theory can then generate the observed (almost) scale 
invariant primordial power spectrum. 



3. SU{2) vector model: equations of motion for the background and for 
linear perturbations 

Let us consider some models where inflation is driven by a scalar field in the presence 
of an 577(2) vector multiplet [97, 98]. A fairly general Lagrangian can be the following 



S = J d 4 x^f^g 



m 2 pR 



E 

a=l,2,3 



rpa rpa 
r liu r a/3 



M 2 
Z a=l,2,3 



(12) 



where is the Lagrangian of the scalar field and = <9 M 77" — d u B® + g c e abc B b ^B^ 
(g c is the SU(2) gauge coupling). Both / and the effective mass M can be viewed as 
generic functions of time. The fields B^ are comoving and related to the physical fields 
by A ^ = (77q, B^/a). The free field operators can be Fourier expanded in their creation 
and annihilation operators 
d 3 q 



5A^x, V ) 



(27T) 



E 



'e*(q)af5Al( q ,r)) + e*\-q) (a^) J 5A?(q,r))], (13) 

X=L,R,long 

where the polarization index A runs over left (7), right (R) and longitudinal (long) 
modes and 



a, A 



(14) 



Here n the conformal time (drj = dt/a(t)). Once the functional forms of / and M have 
been specified, the equations of motion for the vector bosons can be written. For the 
most part, the calculations are quite general in this respect. In fact, the expression of all 
correlation functions, prior to explicitating the wavefunction for the gauge bosons, apply 
to any SU(2) theory with an action as in (12), both for what we will call the "Abelian" 
and for the "non- Abelian" contributions. In particular, the structure of the interaction 
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Hamiltonian is independent of the functional dependence of / and M and determines 
the general form of and the anisotropy coefficients appearing in the final "non-Abelian" 
expressions (see Sec. 4). When it comes to explicitate the wavefunctions, a choice that 
can help keeping the result as easy to generalize as possible is the following 



5B 1 



Js/2(x) +iJ-3/2{x) , (15) 



2y/k 

for the transverse mode and 

5B n = n(x)8B T , (16) 

for the longitudinal mode {n is a unknown function of x = —kt]) [97, 98]. Let us see 
why. As previously stated, for / ~ a a and with a = 0, 1, —2, it is possible to verify 
that the (physical) transverse mode behaves exactly like a light scalar field in a de 
Sitter background. Considering the solution (15) then takes into account at least these 
special cases. As to the longitudinal mode, a parametrization was adopted as in (16) 
in order to keep the analysis more general and given that, because of the instability 
issues, introducing this degree of freedom into the theory requires special attention. 
We are going to keep the longitudinal mode "alive" in the calculations we present, by 
considering a nonzero function n(x), and focus on the simplest case of / = 1. This case 
is known to be affected by quantum instabilities in the longitudinal mode, anyway we 
choose / = 1 for the sake of simplicity in our presentation. The results can be easily 
generalized to gauge invariant models (please refer to [98] for a sample generalization of 
the calculations to massless / ~ a^~ 2 ^ models). 



4. Correlation functions of curvature fluctuations: analytic expressions 

We are now ready to review the computation of the power spectrum, bispectrum and 
trispectrum for the curvature fluctuations ( generated during inflation 

= (27r)V*>(£ + k 2 )P c (k), (17) 

(CkSk&s) = (2*) a 6®fc + k 2 + k 3 )B c (k u k 2 , k 3 ) (18) 

(CfcCfeCSaCfc) = (27r) 3 5 (3) (£i + h + k 3 + h)T c (h, 4 k 3 , h). (19) 

Notice that, on the right-hand side of (17) through (19), we indicated a dependence 
from the direction of the wavevectors; in models of inflation where isotropy is preserved, 
the power spectrum and the bispectrum only depend on the moduli of the wave vectors. 
This will not be the case for the SU(2) model. 
The <5N formula [123, 124, 125, 126] will be employed 

C(x, t) = N(x, t*,t)- N(t*,t) = 5N(x, t), (20) 

which holds if times t* and t are chosen, respectively, on a flat and on a uniform density 
temporal slices (N is the number of e-foldings of inflation occurring between these two 
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times) X. In the presence of a single scalar field </>, Eq. (20) is further expandable as 

AT(n)(f* f\ 

ax,t) = Y, 1 ' J (^,or (2i) 

where is the partial derivative of the e-folding number w.r.t. on the initial 

hypersurface t*. 

If we apply (21) to the inflaton+S'[/(2)vector model, we have 

C(x , t) = + K5AI + l -N H (50) 2 + l -N%5Al5A b u + A^M" 

+ ^AW<^) 4 + -NSgSAtfAlSAiSAi; + (22) 
where now 

and so on for higher order derivatives. 

Our plan is to show the derivation the correlation functions of £ from the ones of 5(j) 
and after a replacement of the <5N expansion (22) in Eqs. (17) through (19). 
The correlation functions can be evaluated using the Schwinger-Keldysh formula 
[46, 47, 48, 49, 50] 



(n|e(t)|ft) = (o 



T (VfoW^] 0/(t) [ T ^-q\H l{ t> )d t>y^ 1 ^ ? (24) 



where, on the left-hand side, the operator and the vacuum Q are in the interacting 
theory whereas, on the right-hand side, all operators are in the so-called "interaction 
picture", i.e. they can be treated as free fields (the Fouries expansion in Eq. (13) thus 
apply), and |0) is the free theory vacuum. 

When calculating the spectra of (, the perturbative expansions in Eq. (22) and (24) will 
be carried out to only include tree-level contributions, neglecting higher order "loop" 
terms, either classical, i.e. from the <5N series, or of quantum origin, i.e. from the 
Schwinger-Keldysh series. Assuming that the SU{2) coupling g c is "small" and that we 
are dealing with "small" fluctuations in the fields and given the fact that a slow-roll 
regime is being assumed, it turns out that it is indeed safe for the two expansions to be 
truncated at tree-level. 



The correlation functions of ( will then result as the sum of scalar, vector and (scalar 
and vector) mixed contributions. As to the vector part, this will be made up of terms 
that are merely generated by the SN expansion, i.e. they only include the zeroth order 
of the in-in formula (we call these terms "Abelian", being them retrievable in the U(l) 

| We employ a spatial metric gij — a 2 (i)e -2 * (e 7 )y, and at linear level the curvature perturbation 
corresponds to Q = — ^ + HSu, where Su is the fluctuation in the total energy density and "J is a scalar 
metric fluctuation. 
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case), and by ( "non-Abelian" ) terms arising from the Schwinger-Keldysh operator ex- 
pansion beyond zeroth order, i.e. from the gauge fields self-interactions. 



Let us now discuss the level of generality of the results we will present in the next 
sections, w.r.t. the choice of a specific Lagrangian. 

The expression for the Abelian contributions provided in Sees, 4.1 and 4.2.1, apply 
to any SU(2) model of gauge interactions with no direct coupling between scalar and 
vector fields (extra terms would be otherwise needed in Eqs. (36) and (37)). The next 
stage in the Abelian contributions computation would be to explicitate the derivatives 
of the e-foldings number and the wavefunctions of the fields: they both depend on the 
equations of motion of the system, therefore the fixing of a specific model is required at 
this point. 

As to the non-Abelian contributions, the results in Eqs. (51) and (52) are completely 
general except for assuming, again, that no direct vector-scalar field coupling exists. 
The structure of Eqs. (67) and (70) is instead due to the choice of a non-Abelian gauge 
group. The expressions of the anisotropy coefficients I n and L n in Eqs. (67) and (70) 
depend on the specific non-Abelian gauge group (for SU(2) one of the I n is given in 
Eq. (69)). Finally, the specific expressions of the isotropic functions F n (a sample of 
which is shown in Eq. (68)) and G n were derived considering the Lagrangian (12) with 
/ = 1 and the eigenfunctions for the vector bosons provided in Eqs. (15) and (16). 

4-1. The power spectrum 

The power spectrum of £ can be straightforwardly derived at tree-level, using the SN 
expansion (22), from the inflaton and the vector fields power spectra 

p c (k) = P is °(k) [l + g ab (k ■ N a ) (k ■ N b ) + is ab k ■ (N a xN b )]. (25) 

The isotropic part of the previous expression has been factorized in 

P is °{k) = N 2 ^(k) + (N c ■ N d ) Pf, (26) 

where we have defined the following combinations 

P± = (l/2)(P^±Pf ), (27) 

from the power spectra for the right, left and longitudinal polarization modes 

Pg =6 ab 6A a R {k,t*)5A b *(k,f), (28) 
Pf =8 ab 8Al(k,t*)8A^(k,t*), (29) 
P± g = 5 ab 5A? ong (k, t*)5A b : ng (k, t* ) (30) 

The anisotropic parts are weighted by the coefficients 



N a N b (P? b - Pf) 

g ab = \ long +^_ (31) 



+ [N c ■ N d ) Pf 
N a N b P ab 



s ab = 7 \ ^ x , (32) 

N$P^ + (N c ■ N d ) Pf 
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(where a sum is intended over indices c and d but not over a and b). Eq. (26) can also 
be written as 

pcd- 



P lso {k) 



cd ' 



after introducing the parameter 



cd 



(33) 



(34) 



Notice that what when we say "isotropic" , as far as the expression for the power spec- 
trum is concerned, we simply mean "independent" of the direction of the wave vector. 
In this case instead, the vector bosons introduce three preferred spatial directions: the 
r.h.s. of Eq. (25) depends on their orientation w.r.t. the wave vector. 

As expected, the coefficients g ab and s ab that weight the anisotropic part of the power 
spectrum are related to /3d, i.e. to the parameters that quantify how much the expansion 
of the universe is affected by the vector bosons compared to the scalar field. 
Assuming no parity violation in the model, we have s ab = 0; the parameters g ab and /3 a & 
are instead unconstrained. In the U(l) case and for parity conserving theories, Eq. (25) 
reduces to [92] 



P ( (k) = P{ so (k) [l + g (k 



n 



(35) 



where n indicates the preferred spatial direction; also one can check that in this simple 
case, if P + ~ P^ and Pi ong = kP + (k ^ 1), the relation g — (k — l)/3/(l + j3) holds, 
where (3 = (Na/N^) 2 (the anisotropy coefficient g is not to be confused with the SU{2) 
coupling constant g c ). If it is safe to assume \g\ -C 1 (see discussion following Eq. (1) 
and references [69, 70]), a similar upper bound can also be placed on (3. 
In the case where more than one special directions exists, as in the SU(2) model, no 
such analysis on the anisotropy data has been so far carried out, the g a parameters 
cannot then be constrained, unless assuming that the three directions converge into a 
single one; in that case a constraint could be placed on the sum \g\ = | J2 a 9 a \i where 



a = 1,2,3 and PAk) 



TDiso 



(k) 



9" 



k ■ n. 



4-2. Higher- order correlators 

We will present the results for the tree-level contributions to the bispectrum and to the 
trispectrum of £• 

These can be classified in two cathegories, that we indicate as "Abelian" and "non- 
Abelian". The former are intended as terms that merely arise from the 5~N expansion 
and are thus retrievable in the Abelian case; the latter are derived from the linear and 
quadratic expansions (in terms of the gauge bosons interaction Hamiltonian) of the 
Schwinger-Keldysh formula and are therefore peculiar to the non-Abelian case. 
We are going to provide both types of contributions, in preparation for discussing and 
comparing their magnitudes later on in Sec. 5. 
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4-2.1. Abelian contributions 

By plugging the 5N expansion (22) in Eqs. (18) and (19), we have 
1 



B ( (k l7 fa, fa) ^ ^N*N H [P^kjP^fa 



perms.] 
K C P (ki)K d Ah)+ perms. 



N+NZNfr \P (f> (fa)Uf u (fa)+ perms. 



for the bispectrum and 

T c (fci, fc 2 , fc 3 , fc 4 ) D N^T^fa, fa, fa, fa) 

+ iVjiV^ [P0(fci)B (|fci + fc 2 |, Afc, ^) + perms. 
+ MiVW 6 \P%(fa)BJk u fa, \fa + fa\) + perms. 



+ N*N t 



(36) 



P^P^P^lfa + fa\)+ perms. 
[P <f> {fa)P 4> (fa)P (j) (fa)+ perms.] 
+ iVX a A% [P;t(fa + fa)P^fa)P^fa)+ perms. 
+ N»KNl c N; d [P^{fa)Pti{fa)P^\fa + fa\)+ perms 
+ JVjiV^iV^ [P <j> (fa)P cf> (fa)P£(fa)+ perms.' 
+ N^N v h N p £ d [P%(fa)P%(fa)P (t> (fa)+ perms. 
+ N^N^N^N" [P^fa)P4\fa + fa\)Pfi(fa) + perms 
+ N^N^N* + k^P^+perms.^ 

+ N^N^N^N*} [P; p c (fci)P^(^)P^(fci + fa) + perms.' 
+ N£NZN>N$ [P£(fa)P%(fa)P%(fa) + perms. 
for the trispectrum. 

Before we proceed with explicitating these quantities and for the rest of the paper, the 
N£ coefficients will be set to zero. This choice was discussed in Sec. 2 and, more in 
details, in Appendix A of [97]. Summarizing, it is possible to verify that the temporal 
mode Bq = is a solution to the equations of motion for the vector bosons, after 
slightly restricting the parameter space of the theory. The adoption of this kind of 
solution, which is related to the assumption of a slow-roll regime for the vector fields, 
implies that the derivatives of iV w.r.t. the temporal mode can be set to zero. 
Let us now provide some definition for the quantities introduced in (36)-(37): we are 
going to switch from the greek indices fi, v, ... to the latin ones, generally used for 
labelling the three spatial directions, in order to stress that all of the vector quantities 
will be from now on three-dimensional 



(37) 



T™ en (k)Pl b + iT° 3 dd (k)P a _ 



T^ifaP^ 



longi 



(3? 
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where 



7~(£) = ef(k)ef(k) + ef{k)ef(k), (39) 



T°f d (k) =i[et(k)e^(k)-er(k)e* H (k)\, (40) 
T^(k) ee e\(k)e](k). (41) 

The polarization vectors are e L (k) = ^ (cos 6* cos — i sin 0, cos 9 sin + z cos 0, — sin 9), 
e R (k) = e* L (k) and e l (k) = k = (sin 9 cos 0, sin # sin 0, cos #), from which we have 

T^ n {k) = 5 l] -k l k v (42) 

=e^A, (43) 

7^ 9 (£) = hh (44) 

The purely scalar terms in Eqs. (36)-(37) are already known from the literature §. As 
to the mixed (scalar-vector) terms, they can be ignored if one considers a Lagrangian 
where there is no direct coupling between the inflaton and the gauge bosons and where 
slow-roll assumptions are introduced for the fields (see Sec. 4 of [98] for a complete 
discussion on this). Let us then look at the (purely) vector part. Its anisotropy features 
can be stressed by rewriting them as follows 

B c (k 1} k 2: h) D ^NiNlN^ik^ih) = M c k N%Mf (45) 

T c (h, 4 4 h) D N^N^N^N^P^P^P^h + h) 
+ N^N^N^N^P^fyP^P^h) 

= M?L%MJ + MlMJM d k Nf d , (46) 

where 

M c k (k) = KP^ik) = Pl c (k) [8 ik Ni + p ac (k)k k (k ■ N a ) + tq ac (k) (k x N a ) k ] (47) 
= NM{k)Nf e 

= P d J(k)[Nl d ■ N l ef + p df (k) (k ■ m (k ■ N l ef ) + iq^(k)k ■ Ni d x N l ef }. (48) 



In the previous equations, we defined 

pac 



pac pac 

P ac (k) ee -J2H— ^L, (49) 



at.: 



r(k) EE jL, (50) 



id) 



with N a ee (Nl Nl Nl) and N'i d ee (JV£, TV; 
Notice that, as for the power spectrum (25), also in Eqs. (45)-(46) the anisotropic parts 
of the expressions are weighted by coefficients that are proportional either to P_ or to 
(Piong — P+)- When these two quantities are equal to zero, the (Abelian) bispectrum 

§ In single-field slow-roll inflation = H^/2k 3 , where is the Hubble rate evaluated at horizon exit; 
the bispectrum and the trispectrum of the scalar field (B$ and T^) can be found in [9, 10, 127, 11, 36] 
(they were also reported in Eqs. (11) and (12) of [98]). 
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Figure 1. Diagrammatic representations of the tree- level contributions to the vector 
fields bispectrum. 



and trispectrum are therefore isotropized. P_ = in parity conserving theories, like the 
ones we have been describing. According to the parametrization (16) of the longitudinal 



mode, we have Pi 



long 



P+ 



{\n{x)f 



DP+. 



4-2.2. Non-Abelian contributions 

We list the non-Abelian terms for the bispectrum 

B c (h, 4 h) ^ KNZN'B^ih, k 2 , k 3 ) 

and for the trispectrum 

T c (h, k 2 , k 3 , h) D KNlN k c N l d T^{K h, h, h) 



(51) 



+ NlNtN+N^ [P^B^ih, k 2 , k 3 + h) + perms. 
+ KNlN k c N l ™ IP^B^ih + k 2 , k 3 , h) + perms. 



The computation of the vector bosons spectra 



(SAfSA'SAl) = 5 (3) (fci + h + fa)B$ 



aba 



ijk ' 



(6 ASSAYS A%6Af) = 6 {3 \h + k 2 + k 3 + k A )T$ 



abed 



(52) 

(53) 
(54) 



will be reviewed in this section. This requires the expansion of the in-in formula up to 
second order in the interaction Hamiltonian 

nf" 



(efo*))Di<T 



drl (H+Jji) - H- nt {ri) ) ) 



(55) 



+ K -^-(T[e / d V ' (HUv ) ~ H- t ( V ')) / drl' (Hf nt ( V ") - H^")) ]). 

— o — OO J — CO 

The interaction Hamiltonian needs to be expanded up to fourth order in the fields 
fluctuations, i.e. H int 



Hint + H int, where 



H 



(3) 



- g c e abc g lk 9 jl (d^B") 5B b k 5B? + g 2 c e eab e eed g ik g1 l B?8B b j 8Bl6Bt (56) 



H^ t = g 2 c e eah e ecd g^g kl 8B^B h k 5B%Bt 



(57) 



Hnt — y c ° c y y u±j i uj - , k uj - , j uj -'i ■ 

To tree-level, the relevant diagrams are pictured in Figs. 1 and 2. By looking at Eqs. (56) 
and (57), we can see that there is a bispectrum diagram that is lower in terms of power 
of the SU(2) coupling (~ g c ) compared to the trispectrum (~ as a matter of fact, 
for symmetry reasons that we are going to discuss later in this section, g\ interaction 
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Figure 2. Diagrammatic representations of the tree- level contributions to the vector 
fields trispectrum: vector-exchange (on the left) and contact-interaction (on the right) 
diagrams. 



terms are needed to provide a non-zero contributions to the bispectrum. 



The propagators for "plus" and "minus" fields are 

5Bt' + {rf)5B^ + W) = ^(V,V"MV - v") + fiJfoVW?" ~ v), (58) 

5Bt' + ( V ; )6B b f(r ] ") = n^( V \ V "), (59) 

SBrW^B^ivl^^i^ri"), (60) 

6B?-(ti')6B}-(ti") = m^ri')Q{ri - V ") + fl%(v\v"Mv" ~ v), (61) 
or 

nj(fc) = T™ en (k)Pl b + iT^{k)P^ + Tl° n9 (k)P% (62) 

n#(£) = T™ en (k)Pl h + iT° 3 dd (k)P% + Tl° n9 (k)P% (63) 



in Fourier space. In the previous equations we set P± b = (1/2) (P$ 6 ± P£ 6 ), P^f = 
5 ab 5Bf(k,r]*)5B*^ b (k,r]) and Pf = (Pf )* (similar definitions apply for Pf and P^ ng ). 

We are now ready to show the computation of the following contributions to the 
bispectrum and trispectrum of ( 

(CfiCsC*) 3 iVXiV c fe (M«(^)54(fe)M c fe (^)>, (64) 

(CfeCfeCfcCs.) 3 ^^(^(^MjC^MgC^H^)). (65) 



Eq. (64) becomes 



^(3)^ + ^ + ^) 



a 3 (?7* 

2 _ea!a_efec 



B? 



drja^flUki^ih^ih) 



+ perms. + c.c. 



(66) 



Even before performing the time integration, one realizes that, because of the 
antisymmetric properties of the Levi-Civita tensor, the ~ g c contribution on the r.h.s. 
of Eq. (66) is equal to zero once the sum over all the possible permutations has been 
performed. The vector bosons bispectrum is therefore proportional to g 2 c . The final 
result from (66) has the following form 
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{CkSkSks) (2vr) 3 5( 3 )(£ 1 + 4 + 4)^ (67) 

n 

where F n are isotropic functions of time and of the moduli of the wave vectors 
(i = 1,2,3) and I n are anisotropic coefficients. The sum in the previous equation 
is taken over all possible combinations of products of three polarization indices, i.e. 
n G (EEE,EEl,ElE,...,lll), where E stands for "even", I for "longitudinal". The 
complete expressions for the terms appearing in the sum are quite lengthy (see Sec. 4.2 
of [97]). As an example, we report one of these terms 

Fin = ~ n 6 (a;*) 2 2 2 [A EEE + (B EEE cosx* + C EEE sinx*) EiX*] (68) 
24k b ki k 2 k%x* z 

hi = e aa ' b 'e ac ' e [( (h ■ N a ') (k 3 ■ N b ') (k 2 ■ N c ') (k y ■ k 2 ) (k ■ A e ) 

- (k 3 ■ N a ) (hi ■ N b ') (h ■ N c ') (h ■ k 2 ) (k ■ A e ) ) + (1 <-»• 3) + (2 3)] (69) 

where A EEE , B EEE and C EEE are functions of x* and of the momenta fc, = \ki\ (they 
are all reported in Appendix C of [97]), Ei is the exponential-integral function and % ■<->• j 
means "exchange fcj with kj" . As we will discuss in more details in Sec. 6, one of the 
more interesting features is that the bispectrum and the trispectrum turn out to have an 
amplitude that is modulated by the preferred directions that break statistical isotropy. 

Let us now move to the trispectrum. Again, we count two different kinds of 
contributions, the first from ~ g c and the second from ~ g 2 c interaction terms, 
respectively in H-^j. and H^ t . The former produce vector-exchange diagrams, the latter 
are represented by contact-interaction diagrams (see Fig. 2). Their analytic expressions 
are different, but they both have a structure similar to (67) 

(dntetft) D (27r) 3 ^ (3) (£i + h + h + h)g 2 c Hl (70) 
x G n (ki, kyii &i4> V )L n (ki ■ kj, Ai ■ Aj, k \ ■ Aj) 

n 

where we define kf 2 = \ki + k 2 \ and kf 4 = \ki + k^\ (see Sees. 5.2.1 and 5.2.2 of [98] for 
the explicit expressions of the functions G n and L n ). 

5. Amplitude of non-Gaussianity: f^L an d t^l 

In this review we use the following definitions for the non-Gaussianity amplitudes 

5 JNL P is °(ki)P iso (k 2 ) + perms. 1 ' 

TNL P iso {ki)P iso {k 2 )P iso {k iA ) + 23 perms. 1 ' 

The choice of normalizing the bispectrum and the trispectrum by the isotropic part of 
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Table 1. Order of magnitude of Jnl in different scenarios. 



NL 



J NL 



fNA 
J NL 



general case 



(1+0)2 ~Nf 



f3 N AA 
(1+13)2 -jvj- 



(1+/9) 



5 5 



v. inflation 



I mpH J 



v.curvaton 



1 + 



1 + 



Am 2 p 

A tot 



1 + 



A 2 r 



7 I AOf 

tot 



the power spectrum, instead of using its complete expression P^, is motivated by the fact 
that the latter would only introduce a correction to the previous equations proportional 
to the anisotropy parameter g, which is a small quantity. 

The parameters f^L and t/vl receive contributions both from scalar ("s") and from 



vector 



fields 

Jnl 
tnl 



As) 
Jnl 

1 NL 



Av) 
Jn Li 

(«) 



r 



A/X- 



The latter can again be distinguished into Abelian (A) and non-Abelian (NA) 



Av) _ f(A) 
JNL 



_(«) 
NL 



Inl 
AM) 
'NL 



, f(^) 

, _(>4a) , -(iVAO (iVA 2 ) 
"r ' ATL "T ' ATL "r 'NL 



The contribution comes from Eq. (45), f^L^ from (67), 
finally from (70) and t^ l A2 ^ from the last line of (52). 



(73) 
(74) 

(75) 
(76) 



and ri A L 2) from (46), 



Notice that, in order to keep the vector contributions manageable and simple in their 
structure, all gauge and vector indices have been purposely neglected at this point and 
so the angular functions appearing in the anisotropy coefficients have been left out of 
the final amplitude results. This is acceptable considering that these functions will in 
general introduce numerical corrections of order one. Nevertheless, it is important to 
keep in mind that the amplitudes also depend on the angular parameters of the theory. 
We will now focus on the dependence of /nl and tnl from the non- angular parameters 
of the theory and quickly draw a comparison among the different contributions listed in 
Eqs. (73) through (76). 



The expression of the number of e-foldings depends on the specific model and, in 
particular, on the mechanism of production of the fluctuations. Two possibilities have 
been described in Sec. 2. For vector inflation we have 

A? Sabft 



2m 2 P ' 



2m 2 P 



(77) 



(see Appendix B of [97] for their derivation). In the vector curvaton model the same 
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Table 2. Order of magnitude of the vector contributions to tnl 
in different scenarios. 



general case 



JVA, 

T ^ 

NL 

1n 3 f 2 ^l I 



NA 2 
NL 



10" 



U) (^) »^A, 



3/2 3/2 2 



T^ 1 
NL 

I™P J, AA 



r A 2 

' JVZ, 

3 3/2 £ 3/2 



.inflation 



same as above 



(i+/3) J 







v.curvaton 



same as above 



10 



(1+/3) J 



r,3 3 / 2 e 3 / 2 /m P \3 
(l + ,3) 3 I A I 



Table 3. Order of magnitude of the ratios f^ L / '/jvx in different scenarios. 





S'nl/ $NL 


fNA / f s 
J NL 1 J NL 


general case 


aN AA 
P N4.4, 


02 2 fmf _^t_ 

P 9c\ H ) n^ 


v. inflation 





P 2 9 2 a (m P \2 
e« V H ) 


v.curvaton 




P 2 sl ( A\ 2 



Table 4. Order of magnitude of the ratios t'n L /t^ l in different scenarios. 





T NAl lT S 

'nl I'nl 


'nl i'nl 


'JVL' 'JVL 


T^ 2 /T 3 

'jvl/ ' nl 


general case 


10 3 /3 2 9c 2 (¥) 2 


10- 5 ^/2 e l/2 g 2 (A) m 2 pNAA 


/3^P^AA 


S/2 e l/2 m 3 pNAAA 


v. inflation 


same as above 


lO -502 e l/2 g 2 (A) (m^ 







v.curvaton 


same as above 


10-*rfPW£ (£) (^) (^i) 2 




r/3 3/2 e l/2 (H^) 3 



quantities become [92, 97] 

AT = ir^L. N'l = ir-*i£- (78) 

" 3 £ t |A<f "* 3 £ C |A=| 2 ^ ' 

Neglecting tensor and gauge indices, the expressions above can be simplified as ~ 
A/m 2 P and Naa — I/Trip in vector inflation, Na ~ r/A and Naa — r/A 2 in the vector 
curvaton model. Also we have Naaa — in vector inflation and Naa a — r/A 3 in vector 
curvaton. 

We are now ready to provide the final expressions for the amplitudes: in Table 1 we list 
all the contributions to /nl, Table 2 includes the vector contributions to t^l, the scalar 
contributions being given by 

r W = "± 3 + $ (79) 

In the expressions appearing in the tables, numerical coefficients of order one have not 
been reported. Also, m is by definition equal to mp in vector inflation and to Aj \fr in 
the vector curvaton model; N^ ~ (mp v /e^)~ 1 and N^ ~ mp 2 , with = (<fi 2 ) / (2mpH 2 ) . 

The quantities involved in the amplitude expressions are g, (3, r, e^, g c , m P /H, A/m P 
and A/H. We already know that g and are to be considered smaller than one (see 
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Table 5. Order of magnitude of the ratios t^ l / (fpf^) 2 m different scenarios. 



NA 1 
NL 



T NL 



NL 



A 
T NL 



NA) 
NL J 



g 2 /3 2 \m P J 



10 



g 3/2 (1 + /3) , ^ 



(— y 







1n 9 rMl+g) mj, h 2 



10 



r 5 e 3 / 2 (l + /3) ff 3 m P mj 



10 



'in 
. -i / 



10 



g r 3 £ 3/2 (1 + ;3) m ^ g4 



~A^~ ~A^ 



discussion after Eq. (35)). Similarly, as mentioned after Eq. (3), r has to remain small at 
least until inflation ends so as to attain an "almost isotropic" expansion. The slow-roll 
parameter and the SU(2) coupling g c are small respectively to allow the inflaton to 
slowly roll down its potential and for perturbation theory to be valid. The ratio mp/H 
is of order 10 5 (assuming ~ 1CT 1 ). Finally, A/mp and A/H have no stringent bounds. 
A reasonable choice could be to assume that the expectation value of the gauge fields 
is no larger than the Planck mass, i.e. A/mp < 1. As to the A/H ratio, different possi- 
bilities are allowed, including the one where it is of order one (see Sec. 6 and Appendix 
A of [97] for a discussion on this). 



Let us now compare the different amplitude contributions. The ratios between scalar 
and vector contributions are shown in Table 3 for the bispectrum and Table 4 for the 
trispectrum. We can observe that the dominance of a given contribution w.r.t. another 
one very much depends on the selected region of parameter space. It turns out that it 
is allowed for the vector contributions to be larger than the scalar ones and also for the 
non-Abelian contributions to be larger than the Abelian ones. This is discussed more 
in details in Sec. 6 of [97]. An interesting point is, for instance, the following: ignoring 
tensor and gauge indices, the ratio g c A/H, that appears in many of the Tables entries, 
is a quantity smaller than one; if we consider the different configurations identified by 
gauge and vector indices, we realize that this is not always true, in fact the value of this 
ratio can be 3> 1 in some configurations. 

Finally, it is interesting to compare bispectrum and trispectrum amplitudes (see table 5). 
Again, it is allowed for the ratios appearing in Table 5 to be either large or small, 
depending on the specific location within the parameter space of the theory. For 
instance, the combination of a small bispectrum with a large trispectrum is permitted. 
The latter is an interesting possibility: if the bispectrum was observably small, we 
could still hope the information about non-Gaussianity to be accessible thanks to the 
trispectrum. 

Another interesting feature of this model is that the bispectrum and the trispectrum 
depend on the same set of quantities. If these correlation functions were independently 
known, that information could then be used to test the theory and place some bounds 
on its parameters. 
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6. Shape of non-Gaussianity and statistical anisotropy features 

Studying the shape of non-Gaussianity means understanding the features of momentum 
dependence of the bispectrum and higher order correlators. If they also depend on 
variables other than momenta, it is important to determine how these other variables 
affect the profiles for any given momentum set-up. This is the far as the 

bispectrum and the trispectrum of the gauge fields are concerned, given the fact that 
they are functions, besides of momenta, also of a large set of angular variables (see 
Eqs. (67) and (70)). 

6.1. Momentum dependence of the bispectrum and trispectrum profiles 

We show the study of the momentum dependence of the F n and G n functions in Eqs. (67) 
and (70) first and then analyze the angular variables dependence of the spectra, once 
the momenta have been fixed in a given configuration. A natural choice would be to 
consider the configuration where the correlators are maximized. 

The maxima can be easily determined for the bispectrum by plotting the isotropic 
functions F n and G n in terms of two of their momenta. These plots are provided in 
Fig. 3, where the variables are x 2 = k 2 jk\ and x$ = k^/ki. Each one of the plots 
corresponds to a single isotropic functions of the sum in Eq. (67). It is apparent that 
the maxima are mostly located in the in the so-called local region, i.e. for k\ ~ k 2 ^> k%; 
three out of the eight graphs do not have their peaks in this configuration but, at the 
same time, they show negligible amplitudes compared to the "local" peaked graphs. 
The situation is much more complex for the trispectrum, being the number of momentum 
variables larger than three (fci, k 2 , k^, k^, fcf 2 and k^). The momentum dependence 
of the isotropic functions can be studied by selecting different configurations for the 
tetrahedron made up by the four momentum vectors, in such a way as to narrow 
the number of independent momentum variables down to two. A list of possible 
configurations was presented in [41]. We consider two of them, the "equilateral" and 
the "specialized planar" . 

In the equilateral configuration the four sides of the tetrahedron have the same length 
(ki = k 2 = k% = k^), therefore x = k 12 /ki and y = k^Jk\ can be chosen as variables for 
the plots. The plots of the isotropic functions of contact interaction and vector exchange 
contributions are provided in Fig. 4. The former (c.i.) shows a constant behaviour in 
this configuration, being independent of k 12 and kf A . The latter (v.e.(I), v.e.(II) and 
v.e.(III)) diverge as k^ 3 [i = 1,2,3 respectively for the three plots) in the limit of a 
flat tetrahedron, i.e. (k^/ki) 0. 

In the specialized planar configuration, the tetrahedron is flattened and, in addition to 
that, three of the six momentum variables are set equal to one another [k\ = k% = k^); 
this leaves two independent variables, which can be x = k 2 jk\ and y = k 3 /ki. There is 
a double degeneracy in this configuration, due to the fact that the quadrangle can have 
internal angles larger than or smaller /equal to tt, as we can see from the plus and minus 
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1.0 "In 




Figure 3. Plot of r n = Q(x 2 — X3)@(x3 — 1 + X2)x\x\R n (x2, 23), where 
we define R n = kfFn. The Heaviside step functions O help restricting 
the plot domain to the region (x2,Xs) that is allowed for the triangle 
k% + /c2 + ^3 = (in particular, we set X3 < X2) • We also set x* = 1. 
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Figure 4. Plots of the isotropic functions appearing in the vector fields trispectrum 
(from Eq. (70)): c.i. is the contribution from contact-interaction diagrams, v.e.(I), 
v.e.(II) and v.e.(III) are the contributions from the vector-exchange diagrams. The 
equilateral configuration has been considered in this figure. 



signs in the expressions for k± 2 and k^ 3 [41] 

k f = ^l + ^±f V^-* 2 )(4-* 2 ), (80) 

k f = f 2 + l/ 2 -^Tf y/(4-a?)(4-V*j. (81) 

The two cases are plotted in Figs. 5 and 6. Notice that divergences generally occur as 
x, y — > 0, as x — > y and (x, y) — > (2, 2). 

6.2. Features and level of anisotropy 

Statistical homogeneity and isotropy are considered characterizing features of the CMB 
fluctuations distribution, if one ignores the issues raised by the "anomalous" detections 
we presented in the introduction. 

Homogeneity of the correlation functions equates translational invariance and hence 
total momentum conservation, as enforced by the delta functions appearing on the left- 
hand sides of Eqs. (17) through (19). This invariance property can then be pictured as 
the three momentum vectors forming a closed triangle for the bispectrum and the four 
momenta arranged in a tetrahedron for the trispectrum (see Fig. 7). 
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Figure 5. Plots of the contact interaction and of the vector-exchange contributions 
in the specialized planar configuration (plus sign). 



Statistical isotropy corresponds to invariance w.r.t. rotations in space of the momentum 
(for the power spectrum) and of the triangle or tetrahedron made up by the momenta, 
respectively for the bispectrum and the trispectrum. This symmetry can be broken, 
as it for example happens in the £77(2) case, by assuming the existence of preferred 
spatial directions in the early universe that might be revealed in the CMB observations. 
When this happens, the correlation functions are expected to be sensitive to the spatial 
orientation of the wave number or of the momenta triangles and tetrahedrons w.r.t. 
these special directions. Analitically, the bispectrum and the trispectrum will depend 
on the angles among the vector bosons and the wave vectors (besides the angles among 
the gauge bosons themselves), as shown in the coefficients appearing in Eqs. (67) and 
(70). This implies that both the amplitude and the shape of bispectrum and trispectrum 
will be affected by these mutual spatial orientations. The modulation of the shapes by 
the directions that break statistical anisotropy was discussed with some examples both 
for the bispectrum and the trispectrum in our previous papers [97, 98]. These examples 
are here reported in Figs. 9 and 10. In Fig. 9 we show the plot of the vector contribution 
to the bispectrum of £, properly normalized in the configuration 

7^ = 7^(0,0,1) (82) 
Ni = N 2 = N^(sm 9 cos </>, sin 6 sin 0, cos 9), (83) 
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where, the (x, y, z) coordinate frame is chosen to be k 3 = x and k\ = k 2 = z and 5 is 
the angle between Ni j2 and k 3 . 

In Fig. 10 we provide a similar plot for the trispectrum, but in a different configuration 

# 2 -^ = 0(1 = 1, ...4) 

N 1 -k 1 = cos 5, N l -k 2 = Q 

N 3 -k 2 = cos 9, N 3 -h = 0. (84) 

In both examples, it is assumed for simplicity that the N a have the same magnitude 
for all a = 1, 2, 3. 

Another comment should be added concerning statistical anisotropy in the model. 
Notice that both the bispectrum and the trispectrum can be written as the sum of 
a purely isotropic and an anisotropic parts. The orders of magnitude of these two parts 
can, for instance, be read from Table 2 for the trispectrum: each one among t^ 2 , 
and Tff% provide the order of magnitude of the level of both their isotropic and 
anisotropic contributions, which are therefore comparable; r^ 1 instead quantifies a 
purely anisotropic contribution which, as discussed in Sec. 5, can be comparable to 
the other three parts, if not the dominant one. A similar discussion applies to the 
bispectrum (see f$ L and f$L i n Table 1). We can then conclude that for the three and 
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Figure 7. Representation of momentum conservation for the bispectrum (the three 
momenta form a closed triangle) and for the trispectrum (the momenta form a 
tetrahedron). 




Figure 8. Plot of f(6,6) = [{B c (9, 5, x* , x 2 , x 3 )xlxlkf )/{glH 2 m 2 N%)} evaluated at 
(x* — 1, X2 = 0.9, X3 = 0.1) in a sample angular configuration. Sec Appendix D of [97] 
for its complete analytic expression. 




Figure 9. Plot of the anisotropic part of the trispectrum from the contribution due 
to vector- exchange diagrams in a sample angular configuration. See Sec. 8 of [98] for 
its analytic expression. 

for the four point function, there is room in the parameter space of the theory for the 
anisotropic contributions to be as large as, or even larger than, the isotropic ones. 
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7. Conclusions 



Motivated by the interest in models that combine non-Gaussianity and statistical 
anisotropy predictions for the CMB fluctuations, we have considered models of inflation 
where primordial vector fields effectively participate in the production of the curvature 
perturbations (. More specifically, we have reviewed the computation of the correlation 
functions up to fourth order, considering an £77(2) vector multiplet. The 5N formalism 
was employed to express ( in terms of the quantum fluctuations of all the primordial 
fields. The Schwinger-Keldysh formula was also used in evaluating the correlators. 

The correlation functions result as the sum of scalar and vector contributions. The latter 
are of two kinds, "Abelian" (i.e. arising from the zeroth order terms in the Schwinger- 
Keldysh expansion) and "non- Abelian" (i.e. originating from the self-interactions of the 
vector fields). The bispectrum and the trispectrum final results are presented as a sum 
of products of isotropic functions of the momenta, F n and G n in Eqs. (67) and (70), 
multiplied by anisotropy coefficients, I n and L n in (67) and (70), which depend on the 
angles between the (gauge and wave) vectors. 

The amplitude of non-Gaussianity has been presented through the parameters /jvl and 
tnl] in particular we have show the dependence of these functions from the non-angular 
parameters of the theory. We have provided the comparisons among the different (scalar 
versus vector, Abelian versus non- Abelian) contributions to /nl and t^l, noticing that 
any one of them can be the dominant contribution depending on the selected region 
of parameter space. In particular, we have stressed how the anisotropic contributions 
to the bispectrum and the trispectrum can overcome the isotropic parts. An interest- 
ing feature of these models is that the bispectrum and the trispectrum depend on the 
same set of parameters and their amplitudes are therefore strictly related to one another. 

We have presented the shapes of both the bispectrum and the trispectrum. The isotropic 
functions appearing in their final expressions had been analyzed separately from their 
anisotropy coefficients. The bispectrum isotropic functions had been found to preferably 
show a local shape. The trispectrum ones had been plotted selecting equilateral and 
specialized planar configurations. The full expressions (i.e. complete of their anisotropy 
coefficients) of bispectrum and trispectrum have been presented in specific momenta 
configuration, in order to provide a hint of the modulation of shapes and amplitudes 
operated by anisotropy. 

We have reviewed old and recent vector field models, indicating both their limits and 
achievements. We would like to stress that, in our view, the most promising features of 
these models consists in the possibility of providing both non-Gaussianity and statistical 
anisotropy predictions that are related to one another because of the fact that they 
share the same underlying theory. This might, at some point in the future, become a 
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great advantage: measurements of non-Gaussianity could be used to constrain statistical 
anisotropy or viceversa. 
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